Astronomy & Astrophysics manuscript no. late 'flares 


©ESO 2010 


October 5, 2010 





Gamma-Ray Burst long lasting X-ray flaring activity 

M.G. Bernardini 1 - 2 , R. Margutti 1 ' 3 , G. Chincarini 1 ' 3 , C. Guidorzi 4 , and J. Mao 1 ' 5 



o 

(N 
-i— > 
O 

O 

X 

43 
6 

(N 
> 

m 
oo 
m 

o 
o 



X 



1 INAF - Osservatorio Astronomico di Brera, via Bianchi 46, 1-23807 Merate (LC), Italy 

2 ICRANet, p.le della Repubblica 10, 1-65 100 Pescara, Italy 

3 University of Milano Bicocca, Physics Dept., Pzza della Scienza 3, Milano 20126, Italy 

4 University of Ferrara, Physics Dept., via Saragat 1, 1-44122 Ferrara, Italy 

5 Yunnan Observatory & Key Laboratory for the Structure and Evolution of Celestial Objects, Chinese Academy of Sciences, 
Kunming, Yunnan Province, 650011, China 

ABSTRACT 

Context. One of the most intriguing features revealed by the Swift satellite is the existence of flares superimposed to the Gamma-Ray 
Burst (GRB) X-ray light curves. The vast majority of flares occurs before 1000 s, but some of them can be found up to 10 6 s after the 
main event. 

Aims. In this paper we shed light on late time (i.e. with peak time t p \ > 1000 s) flaring activity. We address the morphology and 

energetic of flares in the window ~ 10 3 - 10 6 s to put constraints on the temporal evolution of the flare properties and to identify possible 

differences in the mechanism producing the early and late time flaring emission, if any. This requires the complete understanding of 

the observational biases affecting the detection of X-ray flares superimposed on a fading continuum at t > 1000 s. 

Methods. We consider all the Swift GRBs that exhibit late time flares. Our sample consists of 36 flares, 14 with redshift measurements. 

We inherit the strategy of data analysis from Chincarini et al. (2010) in order to make a direct comparison with the early time flare 

properties. 

Results. The morphology of the flare light curve is the same for both early time and late time flares, while they differ energetically. 
The width of late time flares increases with time similarly to the early time flares. Simulations confirmed that the increase of the width 
with time is not due to the decaying statistics, at least up to 10 4 s. The energy output of late time flares is one order of magnitude 
lower than the early time flare one, being ~ l%E pmmp ,. The evolution of the peak luminosity as well as the distribution of the peak 
flux-to-continuum ratio for late time flares indicate that the flaring emission is decoupled from the underlying continuum, differently 
from early time flares/steep decay. A sizable fraction of late time flares are compatible with afterglow variability. 
Conclusions. The internal shock origin seems the most promising explanation for flares. However, some differences that emerge 
between late and early time flares suggest that there could be no unique explanation about the nature of late time flares. 

Key words, gamma-ray: bursts - radiation mechanism: non-thermal - X-rays 



1. Introduction 

One of the most intriguing a nd unexpected feature s revealed by 
the X-Ray Tel escope (X RT, Burro ws et al.ll2005al) on board the 
Swift satellite (iGehrels et alj 2004 ) is the existence of flares su- 
perimposed to t he Gamma-Ray Burst (GRB) X -ray light curves 
(Bur rows et alj l200 5b; Fal cone et alj l2.006h IChincarini et al.l 
2007; Falcone et al.ll2007l). The vast majority of flares occurs 
before 1000 s (Chincarini et al. 2007), but some of them can be 
found up to 10 6 s after the main event. 

Recent analyses of the flare tempor al and spectral prope rties 
( IChincarini et alj|2010l hereafter C10: iMargutti et alj|2010ch of 
a large sample of early time (i.e. with peak time t p k < 1000 s) 
flares and of a subsample of bright flares revealed close simi- 
larities between them and the prompt e mission pulses, pointing 
to an internal origin of their emission (Chincar ini et al ] 120071 
120101) . Therefore, the central engine itself should remain active 
and variable for long time. Alternatively, flaring emission can 
be powered by dela yed magnetic d issipation during the deceler- 
ation of the ejecta dGian nios 2006). Despite many efforts from 
both the observational and the theoretical point of view, a clear 
explanation of the X-ray flares is still missing (C10). 

In this paper we shed light on late time (i.e. t p k > 1000 s) flar- 
ing activity observed by XRT in the 0.3 - 10 keV energy band. Its 
existence poses constraints on the theoretical models since late 



time flares require huge releases of energy (~ 10 50 erg) up to 1 
month after the main event. A previous work on late time X-ray 
flares has been presented bv lCurran et al.l (|2008). These authors 
analysed a sample of 7 GRBs that exhibit flares after 10 4 s by 
fitting them with a gaussian profile plus power-law underlying 
continuum. The purpose was to compare the temporal and spec- 
tral properties of the underlying continuum and of the tempo- 
ral and flux amplitude va riability of these flares with the results 
found in lChincarini et al.l d2007l) and with the prescription of the 
internal and the external shock models. They concluded that late 
time flares are not different from the early time ones. However, 
due to the small number of flares in their sample, this statement 
needed further investigation. 

The present work addresses the morphology of flares in the 
temporal window ~ 10 3 - 10 6 s of a larger sample than in 
ICurran et al.1 d2008). The good fraction of flares of the sample 
with redshift measurement allows us also to characterise the en- 
ergetic properties at very late times. We inherit the C10 strategy 
of data analysis to make a direct comparison with the results ob- 
tained for early time flares and put stringent constraints on the 
temporal evolution of the flare properties. The aim is to identify 
possible differences in the mechanism producing the observed 
early and late time flare emission, if any. This requires the com- 
plete understanding of the observational biases affecting the de- 
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tection of X-ray flares superimposed on a fading continuum at 
t > 1000 s. 

In Sect. [2] we present the late time flare sample and we de- 
scribe the fitting procedure. In Sect.[3]we show the main results 
of the present analysis of the late time flares temporal profiles 
and we compare them with those obtained from the CIO sam- 
ple. In Sect. |4] we discuss our findings. In Sect.[5]we summarize 
the results obtained and draw our main conclusions. We adopt 
standard values of the cosmological parameters: H a = 70 km 
s _1 Mpc~', Q.m = 0.27 and Qa = 0.73. Errors are given at 1 cr 
confidence level unless otherwise stated. 



2. Sample and fitting procedure 

We consider all the Swift GRBs observed between April 2005 
and December 2009 that exhibit late time (i.e. t p k ^ 1000 s) 
flares. Among all the X-ray light curves observed by XRT, we 
include in our sample only those that contain a flaring activity 
with a relatively complete structure: rise, peak and decay phase. 
The present data set consists of 22 long GRBs (LGRB) and 1 
short GRB (SGRB), GRB050724, with 36 flares (see Table ??). 
9 GRBs (14 flares) have measured redshiffl These flares and 
their underlying continuum have been analysed in the 0.3 - 10 
keV energy band. For a full description of the data reduction we 
refer to C10. 

The flare properties have been investigated by fitting the 
0. 3 - 10 keV l i ght curve with an empirical function proposed 
bv lNorris etaP 00051 hereafter Norris05): 




Fig. 1. Width vs. peak time for the sample of late time flares 
(blue squares) compared with early time flare sample (red stars). 
Green line: linear (w = 10. + 0.2 t P k) behaviour found in the 
early flare sample. Pink line: best-fitting for the late flares sam- 
ple, w = 10 ( -°- 9±a7) f^ 1±a2) . Gray area: region of the w - t pk 
plane that cannot be explored in our analysis. Inset: distribution 
of the ratio w/t P k of late time flares (blue rectangles) compared 
with the early time flare one (red rectangles). The two distribu- 
tions are centered on a similar value (the median is 0.28 for the 
late time flares and 0.23 for the early time flares) but late time 
flares have a larger standard deviation (cr = 0.33 compared with 
o- = 0.14). 



I{i)=AAe~ 



for t > u 



(1) 



where fi = (ti/ti) 1 ^ 2 and A = e 2f *. The intensity peaks at t p k = 



T p k + fj — (t\T2) 1 + t s . The pulse width is measured between 
the two l/e intensity points, 



w = At 1/e = t d + t r = t 2 (1 +4(i) 



1/2 



The pulse asymmetry is 



td+tr 



(2) 



(3) 



The decay time and the rise time t r are expressed in terms of 
w and k as: 

t dtr = ^w{\±k). (4) 

The choice of the Norris05 profile allows us to perform a direct 
comparison of the properties of late and early time flares. 

We fitted simultaneously the underlying continuum and the 
flares by adopting a multiply broken power-law plus an appro- 
priate number of Norris05 functions. As a measure of the quality 
of the fit we used the^ 2 statistics. The results of the fitting pro- 
cedure are available as online material. The errors on the derived 
quantities are obtained accounting for the entire covariance ma- 
trix of the fitting parameters. The flare fluence S is calculated by 
integrating the corresponding Norris05 function over the entire 
duration. The count -to -flux conversion factor has been com- 
puted fro m a spectrum extracted during the entire d uration of 
the burst (Evan s et alJl2007l 120091: iMargutti et all2010bl) . 



1 For the redshift values 
http://www.mpe.mpg.de/~jcg/grbgen.html 



refer 



to 



3. Analysis and results 

In the following we characterise the late time flare emission by 
investigating the flare morphology and energetics, that the theo- 
retical models most commonly address. The morphology of the 
flare light curves is presented in the observer frame, while the 
energetic properties are investigated in the source rest frame. 



3.1. Morphology 

Flares are erratic emission of huge amount of energy on a rela- 
tively short timescale. The median of the distribution of the ratio 
between the width w and the time t p k at which flares occur is 
< w/t P k >= 0.28, that is similar to the one of the population of 
early time flares (0.23). However, late time flares have a much 
broader distribution (cr = 0.33 compared with <x - 0.14). The 
highest w/t p k occurs for the flare belonging to the only SGRB of 
our sample (GRB050724). 

In Fig. Q] we show the evolution of w with time for late time 
flares. The width depends on time as w oc t l pk l . The increase of the 
flare duration with time observed in this sample as well as in the 
early time flare sample is one of the main differences from the 
prompt emission pulses: while w is almost constant with time for 
the BATSE pulses belonging to the Norris05 sample, it linearly 
increases for the early time flare sample (see Fig. 3 in C10). 

An important issue is if this increasing trend is a real feature 
or it is due to the progressive widening of the light curve tem- 
poral bins: as the time proceeds, the decaying light curve is built 
by integrating the observed counts over larger time intervals in 
order to keep the signal to noise ratio constant. Therefore, an in- 
tense and narrow flare would be spread into a less intense and 
wider one, thu s resulting in a possibl e bias for the correlation. 
Simulations in lChin carini et ail d2007l) revealed indeed that very 
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short duration flares (up to a factor 1CT 2 with respect to the best 
fitting relation shown in Fig. [TJ could be in principle detected at 
10 5 s at 90% confidence level. This result has been obtained by 
generating flares with any possible luminosity, width and tempo- 
ral occurrence. The real situation is slightly different, since flares 
become dimmer as the time proceeds (see Sect. 13.21 1. This mo- 
tivated us to perform simulations that reproduce more in detail 
the observed scenario. 

We simulate for each flare of the sample a set of fake flares 
by adding to the best fitting of the continuum different Norris05 
profiles that have the same peak intensity and peak time of the 
original flare but with different decreasing widths. Then, we re- 
bin each simulated light curve as the corresponding observed 
one, and we fit the results with a new Norris05 profile plus the 
underlying continuum. We define in this way a minimum observ- 
able width for each time interval as the width calculated with the 
best fitting Norris05 parameters of the shortest simulated flare 
occurring in that time interval that we can identify after the re- 
binning procedure. Low significance fluctuations (number of cr 
< 2) have not been classified as flares. The result is shown in Fig. 
Q] the gray area represents the forbidden region for flare detec- 
tion, namely a flare with those width and peak time values would 
be classified as a non significant fluctuation around the underly- 
ing continuum. We found that the w - t p k best fitting relation is 
unbiased up to 10 4 s since the forbidden region is at least 1 order 
of magnitude below th e best fitting relation, consistently with 
IChincarini et al.l (120071) . After ~ 10 4 s the range of detectable 
widths reduces progressively, being only 24% of the observed 
value at 10 6 s. This is due to the progressively lower statistics of 
the GRB afterglows with time, which possibly causes a depar- 
ture from the best fitting w - t p k relation found at early times. 

Another source of bias affecting in some cases the late time 
flare sample is the difficulty to distinguish blended structures. 
In fact the integration of the observed counts over progressively 
larger time intervals makes harder and harder to discriminate if 
the observed excess is due to one single flare or it comes from the 
contribution of many overlapping flares. Therefore in some cases 
the width obtained with our fitting procedure could be overesti- 
mated but only by a factor of a few. This factor will be not so 
different from the average number of flares observed in a tempo- 
ral interval At ~ t of the early time sample, that is ~ 1.5. 

A remarkable property of the flare profile is a self-similar be- 
haviour, being the rise over decay time constant in time. This fea- 
ture is found in the present sample of late time flares as well as in 
the CIO analysis, with the same best fitting relation td ~ 2 t r and 
a distribution of t r /td (see Fig. |2]i and of asymmetry k (median 
value 0.3 1 with standard deviation 0.22) with a similar central 
tendency and dispersion. The K-S test on the two distributions 
of t r /ta for early time and late time flares gives a 60% probabil- 
ity that they belong to the same population. Still, there is no flare 
having t r > td (gray area in Fig.|2]l. 

3.2. Energy and luminosity 

In Fig. [3] Inset the distribution of the logarithm of the isotropic 
energy emitted in the observer frame 0.3 - 10 keV energy band 
of late time flares is shown, compared with that obtained in C10 
for early time flares. Late flares are less energetic, with a median 
value of 1.4 x 10 50 erg, one order of magnitude lower than the 
early time flare sample. As stated in C10, the distribution func- 
tion of the isotropic energy emitted is likely incomplete at the 
low-energy tail. On the other hand it is evident that high ener- 
getic flares are not detected at late times (see Fig. [3]). 




1 2 3 4 5 6 7 



Log [t r /s] 

Fig. 2. Decay time vs. rise time for the sample of late time 
flares (blue squares) compared with the early time flare sam- 
ple (red stars). The best linear fitting for the two samples is 
t d = (-0.1 + 0.6) + (2. 15 +0.04) t r (pink line). Gray area: t d < t r . 
Inset: distribution of the ratio t r /td of late time flares (blue rect- 
angles) compared with the early time flare one (red rectangles). 
The two distributions are centered on a similar value (median 
0.52 with cr = 0.23 for the late time flares and median 0.49 with 
cr = 0.26 for the early time flares). 



53 




Fig. 3. Isotropic energy emitted in the observer 0.3 - 10 keV en- 
ergy band vs. rest frame peak time for the sample of late time 
flares (blue squares) compared with the early time flare sample 
(red stars). The absence of energetic late time flares is manifest. 
Inset: distribution of the logarithm of the isotropic energy emit- 
ted in the 0.3 - 10 keV energy band of late time flares (blue rect- 
angles) compared with the early time flare one (red rectangles). 
The late time flares are clearly less energetic than early time 
ones, with a median for the distribution of the energy 1.4 x 10 50 
erg and standard deviation cr = 3.2 x 10 50 erg. 

One of the most important clue on the understanding of X- 
ray flares is the evolution of the peak luminosity L p k versus 
time. In Fig. [4] we compare the behaviour L p k °c t~f found in 
C10 with our sample: the flares with > 10 4 s differ from 
this relation for more than 2cr. The best fitting of the late time 
flare luminosity-time relation leads to L p k oc t~ l k 2 . This result 
is consiste nt with the average lum inosity behaviour after 1000 
s found in Margut tTeTall (I2010al hereafter M10). In M10 is 
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Log [r pt "7s] 



Fig. 4. Peak luminosity vs. peak time measured in the rest frame 
for the 14 late time flares with measured redshift (blue squares) 
compared with the early time flare sample (red stars). The best 
fitting for the whole sample accounting for the extrinsic scatter 
is Log[L pk ] = (53.2 + 0.3) + (-1.8 + 0.1) Log[t pk ] with cr ex , = 
(0.72 + 0.07) (pink line). The best fitting for early time flares 
accounting for the extrinsic scatter LoglLpk] = (54.8 + 4.0) + 
(-2.7 + 0.5) Log[t pk ] with cr ext = (0.73 + 0.08) (green line) is 
shown for comparison. The thinner lines indicate the 2cr regions 
for the two best fitting relations. 




-2-10 1 2 3 

Log[DF/B] 

Fig. 5. Distribution of the logarithm of the peak flux to contin- 
uum ratio of late time X-ray flares (blue rectangles) compared 
with the early time flare one (red rectangles). The late time flares 
have clearly a smaller AF/F than early time ones, with a median 
for the distribution of the logarithm 0.04 compared with 0.60 for 
the early time flares, and standard deviation cr - 0.75 compared 
with a- = 0.73. 



clearly shown that at this epoch the average luminosity function 
of flares tracks the detectability threshold. Therefore, the unbi- 
assed luminosity-time relation for late flares is likely steeper. 
If we perform a joined fitting of the early and late time flare 
data, we find L pk oc t~l* '. If we account for an extra-variance 
of the sample (D'Agostini 2005), this best-fitting relation has a 
X 2 /dof = 58.1/55 and a p-value of 40%. 

In Ml we reported that the distribution of the peak flux-to- 
continuum ratio AF/F of the C10 early time flare sample can be 
fitted by the superposition of two gaussian profiles. If we include 
the late time flares of the present sample, we obtain that they all 



lie in the "fainter" part of the distribution (see Fig.[5j. A K-S test 
performed on the two distributions gives a 2 x 10~ 5 probability 
that they belong to the same population. 

The smaller ratio at late times can be ascribed to two occur- 
rences: a) the peak flux of flares decreases with time, then at late 
times flares are fainter (see also Fig.[4]l; b) after 1000 s in 80% of 
the light curves of our sample the underlying continuum changes 
its behaviour becoming shallower. To investigate these two ef- 
fects, we calculate how the AF/F of a late time flare would ap- 
pear at early times. The ratio AF/F is substituted for the equiv- 
alent peak luminosity to continuum ratio AL/L. We assume that 
the peak luminosity of flare evolves as r 2J for t < 1000 s and t~ a 
for t > 1000 s, while the continuum evolves as t for t < 1000 
s (M10) and r l for t > 1000 s, this index being the average 
behaviour between 1000 s and 10 4 s of the light curves of our 
sample. We leave a as a free parameter due to the impossibility 
to evaluate quantitatively the luminosity-time relation for late 
time flares. However, it is reasonably between the value of the 
detection threshold found in M10 and the value found for early 
time flares, 1.2 < a < 2.7. We consider a flare with AL/L = 1, 
equal to the median value for the late flare sample distribution, 
at 10 4 s. Extrapolating backwards, we find that, depending on a, 
the corresponding AL/L at 10 2 s would be larger up to 50, within 
the range of values observed for the early time sample. This sim- 
ple calculation reveals that the discrepancy in the AL/L values 
for late time and early time flares depends on the shallower be- 
haviour of the continuum after 1000 s that results in a bias for 
the flare detection (M10). 

4. Discussion 

The existence of flares at 10 6 s after the trigger poses several 
constraints on GRB central engines. Many authors pointed out 
that late time flares must require late ejecti on of shells rather 
than late collisions of prompt emission shells (|Zhang et al.l 2006 ; 
lLazzati & Pernal 120071: iMaxham & Zhangl [2009). This conclu- 
sion, based o n the constraints impo sed by the observed dura- 
tion of flares (lLazzati & Pernall2007l) . implies that the collisions 
producing flares occur in the same conditions as those produc- 
ing the prompt emission pulses, but th e ejection time is com - 
parable to the flare observational time (lLazzati & Pernal 120071) . 
Therefore, the differences observed in flares with respect to the 
prompt emission pulses depend entirely on intrinsic differences 
in th e shells ejected by the cen tral engine and not on their evolu- 
tion dMaxham & Zhana l2009). More important, the central en- 
gine itself should remain active up to 10 6 s. 

Late time flares are less energetic than early flares by at least 
1 order of magnitude (see Fig. [3] Inset). In Fig. [3] it is mani- 
fest that early time flares from C10 span a wide range of ener- 
gies, and a hint of bimodality appears in their energy distribution 
(see C10, Fig. 8, Inset). Late time flares belong mainly to the 
low energy tail of the distribution. This translates, within the in- 
ternal shock model, into less energetic collis ions at late times 
( lLazzati et all 120081: IMaxham & Zhangl 120091) . and that corre- 
sponds to l ess energetic shells e jected at late times from the cen- 
tral engine (lLazzati et al.ll2008l) . 

The peak luminosity of early time flares decreases with time 
as a power-law L pk oc rf 1 , as found in C10 (see Fig. |4j and ob- 
tained also in M10 with a different method. The late time flares 
of the present sample observed at 10 3 < t < 10 4 s are slightly 
brighter with respect to the expectations based on the early time 
flares, but still consistent within the dispersion (see Fig. |4j». This 
is mainly a result of the bias introduced by the underlying con- 
tinuum: at t ~ 1000 s we observe in 90% of the X-ray light 
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curves of our subsample of GRBs that show late time flares and 
with a redshift measure ment the transition from a steep decay to 
a shallow decay phase dEvans et al.ll2.Q09h . and this prevents the 
flares belonging to the faint tail of the distribution from being 
detected. 

A different situation occurs for the three flaresEI observed at 
t > 10 5 s that are well above the extrapolation of the t~ 2t be- 
haviour. To understand the nature of these flares we estimate 
the probability that they are compatible with the distribution 
of the observed early time flares around the extrapolated fit- 
ting L p it oc r 2J with the same extra-variance of the sample 
o~ext = 0.73. We compute the logarithmic average X-ray lumi- 
nosity light curve of the underlying continuum of the GRBs of 
our sample, and we define the minimum luminosity for a flare at 
a given time as the one that can be detected at 2<x over the av- 
erage continuum light curve. At t — 10 5 s the probability that a 
flare belonging to that distribution has a luminosity greater than 
the minimum one is less than 10~ 5 . At t — 10 6 s this probability 
is even much lower, < 10~ 9 . 

Conversely, we could assume that all the flares of the CIO 
and of the present sample follow the same luminosity-time re- 
lation, but it is not the one found in CIO (and, we recall, also 
in M10 with a different method). In this case for the joined 
sample of early time and late time flares it reads: Log[L] = 
53.2 - l.8Log[t], with an extrinsic scatter <j ext = 0.72. This is 
an acceptable fit (x 1 /dof - 58.1/55, p-value= 40%). However, 
in Ml we showed that the early time flare average luminosity 
between 60 s and 400 s is steeper than r l & at 90% confidence 
level. Moreover the continuum after 1000 s introduces a severe 
bias for the detectability of faint late time flares: probably only 
the brightest flares can overshine the continuum at late times 
(M10). Therefore, the best fitting t~ l 8 likely overestimates the 
peak luminosity of flares. 

If we consider the distribution of the peak flux to continuum 
ratio AF/ F for late time and early time flares, we find that the 
two datasets differ significantly, with a K-S probability that they 
belong to the same population of 2 x 10~ 5 . The faintness of the 
AF/F for late time flares can be ascribed to the shallower be- 
haviour of the underlying continuum after 1000 s: with a simple 
calculation we show that the same flare with AF/F ~ 1 at 10 4 s 
would have a AF/F ~ 50 at 10 2 s due to the steeper continuum. 
This under the crucial assumption that the evolution of the peak 
luminosity of late time flares is decoupled from the underlying 
continuum, differently from what has been found for early time 
flares. In M10 we found that the average luminosity of flares at 
t < 400 s decreases with time with a power-law index that is 
linearly correlated to the one of the underlying steep decay. This 
points to a deep connection between the mechanisms that pro- 
duce the two components. If this connection is not established at 
late times, it is natural to expect that late time flares are so rare, 
since only the brightest ones are able to overshine the continuum 
(M10). 

In the framework of accretion models, the long-lasting ac- 
cretion up to late times can be due either to a very small vis- 
cosity parameter a < 0.01, or to a continued mass fall-back 
dKumar et al.1 12008). The first possibility seems to be disfavored 
since it predicts temporal indices for the initial and final de- 
cays of the whole li ght curve less steep than the observed ones 
(Kum ar et all |2008). In the second scenario, in a very short 
timescale a quasi steady-state is established with m acc ~ mft,, 
where m acc and m/j are the accretion rate and the mass fall- 



2 Note however that one of the 3 flares with t pk > 10 5 s belongs to the 
short GRB050724. 



back rate respectively; therefore t he observed lumino sity tracks 
the mass fall-back rate fairly well dKumar et al.l l2008). The fall- 
back radius is rji, > 1000 r s , where r s is the corresponding 
Schwarzschild radius (for a central black hole mass Mbh - 
5M ), and an accretion timescale of ~ 10 4 s is no t unreason- 
able fo r very low accretion rates dPerna et al. 1 12006b . that is the 
case of iKumar et"all (120081) . Fig. 7. 

In this context, the presence of flare s in the light curve is due 
to sudden changes in the accretion rate ( IKumar et alj |2008). The 
duration of each accretion episode, that corresponds t o the dura- 
tion o f the flare, is of the order of the accretion time ( Perna etalJ 
2006) with t acc /t often less than unity (see IKumar et alj 2008, 
Fig. 6-7). Disk instabilities could in principle lead to large varia- 
tions in the accretion rate. A promising possibility is gravitation 
instabilities (we refer to M10 for a discussion on thermal and vis- 
cous instabilities in the context of early time flares): these insta- 
bilities may occur in the outer regions of the disk (r > 100 r s ) if 
the angular momentum is high. This requirement is more likely 
fulfil led in collap sars than in the merging of compact objects 
dPerna et al.l2006l) . Alternatively, the fall-back material can frag- 
ment into bound objects pr odu cing a modulati on in the accre- 
tion rate (Kumar et al. 2008; see Rosswog 2007 in the context of 
the merging of compact objects). If fragments form, their initial 
mass depends strongly on the local properties. However, their 
estimated mass is Mf rag ~ 1 50 erg for fragm entation taking 
place at very large radii/times dPerna et alj 2006. eq. 5), that cor- 
responds to the median energy observed for late time flares (see 
Fig.g] Inset). 

A different process of fragmentation can occur in the col- 
lapse of a rapidly rotating stel lar core, with th e formation of 
more than one compact object (King et al. 2005). These subse- 
quently coalesce under the effect of gravitational radiation re- 
leasing an energy ~ 2% of the main event with a dela y of hun- 
dreds of seconds dDavies et al.ll2002T: iKing et a"Dl2005l) . The re- 
activated burst has a jet axis close to the first one, thus it expands 
in a cleaner environment with lower baryonic contamination an d 
consequently with an higher Lorentz factor (IKing et al.l l2005). 
This can explain why flares are softer than the main event. 

Many GRBs have X-ray light curves without flares. In fact a 
flaring activity is observed in ~ 33% of XRT light curves. The 
conditions for fragmentation in the outer layer of the d isk are un- 
certai n, however high rotation is generally required dPerna et al.l 
2006). Rotation plays an import ant role in shaping the X-ray 
light curve at very early stages dKumar et al.1 12008): Fig. 5 of 
IKumar et al.l (J2008) shows clearly that lowering the rotation rate 
of the central object, the resulting light curve exhibits a longer 
and fainter prompt duration with a steeper decline. However the 
rotation rate cannot be very low (less tha n 8 times lower than 
the model in Fig. 2 of IKumar et al.ll20 08) otherwise the entire 
core directly collapses to a black hole preventing the occurrence 
of the GRB. We noted in M10 that GRBs with multiple flares 
have generally a much flatter steep decay light curve. A flatter 
decay can be produced by a higher rotation rate for these GRB 
central objects, that allows to produce more suitable conditions 
for fragmentation. 

For late time X-ray light curve the angular momentum of the 
fall-back material is crucial to shape the transition between the 
shallow decay phase and the following normal decay. If the an- 
gular momentum is low, this transition results in a sudden drop 
of the light curve. In the opposite case the transition is smooth. 
In our sample of GRBs with late time flares, only 3 cases show a 
final decay steeper than 1.5. This is consistent with the expecta- 
tions of high angular momentum for the fall-back material, and 
to the leading role that it has in favoring the fragmentation. A 
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complete analysis of the light curves without flares needed in or- 
der to confirm this scenario is beyond the scope of the present 
work and will be presented elsewhere. 

Magn etic fields strongly modi fy the dynamics of accretion 
(see e.g. iProga & Begelmanll2003l) . In particular the energy re- 
lease can be repeatedly stopped and restarted by the magnetic 
flux accumulated around the central object, giving origin to the 
obser ved X-ray flares, if the mass supply rate decre ases with 
time (Pro ga & Begelmanl 120031: IProga & Zh ang 2006). In this 
context a flaring episode depends on the strength of the magnetic 
field, namely when the magnetic stress provides support against 
gravity: 2BrBz/4-7t ~ GMS/r 2 , where B is the magnetic field, 
M the mass of the central object, r the radius and 2 the surface 
density. A magnetic flux of the order of (p — nr 2 B{r) ~ 10 29 cm 2 
G is necessary to stop accretion with rate 10~ 3 M s _1 around 
a 3 M central object at r = 300 r^, representative of the late 
time evolution dProga & Zhangl 12006). Simple estimates show 
for the emission episodes that they are less energetic of one or 
two orders of magnitude than the prompt emission, and the du- 
ration of each episode is coupled with the time of occurrence 
(Proga & Zhang 2006), thus reproducing the main observational 
features of flares. 

Despite the similarity with the early time flares and the 
prompt emission pulses points to an internal shock origin of 
late time flares, thus favoring accretion models, there are differ- 
ossibilities as the magnetic reconne ction model 



ent promising possibilities as the magnetic reconne ction model 
dLvut ikov 20061 iGianniosi 120061 lLazar etalJl2009h that do not 
require a long lasting activity of the central engine. In fact in 
this model the emission takes place just before the deceleration 
radius. Within this context, flares are originated by del ayed mag- 
netic dissipations during the deceleration of the ejecta dGianniosI 
2006). For a constant variability (which is indeed the case of 
flares) the isotropic en ergy emitted is p roportional to the energy 
in the forward shock dGianniosI [2006) . This is consistent with 
the fact that late time flares are on average less energetic than 
early time ones. However, the time dependence of the width can 
be recovered only with ad-hoc assumptions on the character- 
istic length scale s of the reconnection regions dGiannio 
lLazar et al.f2009l) . 

The linear evolution of the decay time as a function of the 
rise time is present in the early as well as in the late time flare 
sample, being tj ~ 2 t r (see Fig. [2). This is a result that is valid 
over 5 orders of magnitude, and it is also property common to the 
prompt emission pulses as shown in C10. This is a very strin- 
gent requirement for any theoretical model, since the rise and 
decay times depend eithe r on the dynamic of the shocks within 



the internal shock model ( Kobava shi et al. 19971). or on the mo 



tion of turbulent eddies or mini-jets (Narayan & Kumar! [2009; 
lLazar et al.ll2009h . 

4.1. Flares associated to variability in the GRB afterglow: the 
case of GRB050724 

In the above discussion we focused mainly on an internal origin 
of flares. However, it is of some interest to investigate late time 
flares in the context of variability of the GRB afterg low. A use- 
ful tool in this respect is the plot proposed by lloka et al.l (2005) 
where they trace the regions of allowance for the increase of the 
flux with respect to the continuum as a function of the temporal 
variability for bumps in the afterglow on the basis of kinematic 
arguments. In Fig. |6]we present our sample of late time flares in 
a "Ioka plot", compared with the results presented in Fig. 13 of 
C10. As stated in Sect. 13. li the range of observed w/tpk of the 
two sample is comparable, while as we noted in Fig. [5] the flux 




-1.0 -0.5 
Log [w/fjk] 

Fig. 6. Increase of the flux with respect to the underlying contin- 
uum vs. temporal variability for the sample of late flares (blue 
squares) compared with the early flare sample (red stars). The 
black stars and square correspond to the SGRB early time and 
late time flares, respectively. The solid lines limit the kinemati- 
cally allowed regions for density fluctuations on-axis (blue line) 
and off-axis (red line), for off-axis multiple regions density fluc- 
tuations (green line), for patchy shell model (blac k line) and for 
refreshed shocks (pink line), see lloka et al.l d2005l) for details. 



increase with respect to the continuum covers a different range 
in the case of late time flares. 

All the late time flares of our sample, with the exception of 
that of GRB050724, fulfill the general requirement for internal 
shocks Af/f < 1. However, ~ 14% are compatible with internal 
shock origin only. Differently from the early time flares, ~ 25% 
of late time flares are compatible also with bumps due to on- 
axis density fluctuations. This reveals that, although the internal 
shock origin is the most promising explanation for flares, a siz- 
able fraction of late time flares can have a different origin than 
early time ones. Therefore, there could be no unique explanation 
about the nature of l ate time flares. Similar results were found by 
ICurran et al.l (|2008) that compared the flare-to-continuum flux 
ve rsus variability of the 8 flares of their sample^with the sample 
of C hincarini et al.ld2007l) . They found that one flare (i.e. ~ 12% 
of their sample) can be explained only with internal shocks, and 
all the other flares can be interpreted both by internal and by ex- 
ternal shocks. However, they do not have any flare compatible 
with on-axis density fluctuations. 

The flare of the short GRB050724 is the most challenging. 
As pointed out in Sect. 13. H it has the longest duration among the 
flares of the joined early and late sample, w/ t p k ~ 1 .4. This trans- 
lates in the Ioka plot in being the only case compatible with the 
geometrical constraint for the multicomponent jet model, the re- 
freshed shock model and marginally with the density fluctuation 
of many regions. 

GRB050724 has also been observ ed at the epoch of the flare 
in the optical and radio frequencies (Male sani et al.ll2007l and 
references therein). The simultaneous analysis of X-ray and op- 



3 We recall that the width in the Norris05 profile is measured at ~ 
37% of the maximum. If we consider the corresponding width measured 
on a Gaussian profile, it is w gmlss = 2.83cr, while the Full Width at 
Half Maximum (FWHM) is w FW hm = 2.35tr. Therefore, the position 
on the x axis of the flares in the Ioka plot crucially depends on the 
profile adopt ed and on the definition of the duration. In the work of 
ICurran et all {2008) the duration of flares is the FWHM. 
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tical data led Males ani et al.l (120071) to conclude that they likely 
belong to the same spectral component (with fiox = (0.72 + 
0.04)) that implies a similar temporal behaviour. Therefore, de- 
spite the poor sampling of the optical light curve, they identify 
the presence of a si multaneous flare in the optical energy band 
(Mal esani et alj|2007l) . A rebrightening in the radio light curve 
is also present, but it is delayed w ith respect to the X-ray and 
optical peaks (Mal esani et al.ll2007l) . 

The association of this flaring emission in different energ y 
bands with the afterglow has been tested by iPanaitescul J2006). 
He included the flare in the forward shock emission in an in- 
terstellar medium environment as an episode of energy injection 
and simulated the observed X-ray, optical and radio light curves. 
He found agreement with the observations, how ever using dif- 
ferent assumptions and a smaller dataset than in Mal esani et alJ 
(l2007h . In the same context of the external shock model, the ob- 
served duration of the flare and spectral index can be accommo- 
dated jnjheadiabaticshock propagating in a wind environment 
dLazzafi & Pernall2007l Fig. 2). 

The simultaneous detection of a flare at different wave- 
lengths is not a f eature common to all th e X-ray flares (see e.g. 
CIO). Moreover, Margutti et al. (2010c) found in two cases of 
very bright flares that during the emission episode, the peak en- 
ergy crosses the X-ray energy band with a temporal dependence 
very close to the temporal decay of the observed flux. Although 
this has not been demonstrated to be a c ommon property of all 
flares, the optical-X-ray SED reported bv lMalesani et al.l d2007l) 
seems to exclude that the peak frequency crosses the X-ray band. 

The above reasoning reveals indeed that the flare observed in 
the X-ray light curve of the short GRB050724 can be interpreted 
within different frameworks. Also the early time flares associ- 
ated to SGRBs included in the C10 sample are all compatible 
with afterglow variability (see Fig.[6j>. While consensus has been 
reached on the origin of flares associated to long GRBs, there is 
no overwhelming evidence on the nature of the flares in SGRBs. 
A devoted study is necessary in order to increase our knowledge 
on the SGRB progenitors and their environment. 

5. Summary and conclusions 

We analysed 23 GRB X-ray light curves observed by XRT that 
exhibit late time (i.e. t p % > 1000 s) flares. The present work in- 
herits the C10 strategy of data analysis but extends the temporal 
window of investigation of three orders of magnitude, from 10 3 
s to 10 6 s, with a larger samp le than in previous works on late 
time flares dCurran et al.l l2008h . We found that late time flares 
are similar to early time flares for the following reasons: 

- the width of late time flares increases with time similarly to 
the early time flares. This is a real property of flares, at least 
up to 10 4 s; 

- the self-similar behaviour of the flare profile, being the decay 
proportional to the rise time, is also a property shared by 
both early time and late time flares. This is a more stringent 
requirement than being simply "asymmetrical pulses". 

However, they differ for the following characteristics: 

- the median energy output of late time flares is one order 
of magnitude lower than the early time flare one, being 

- the bias introduced by the shallower underlying continuum 
after 1000 s allows the detection of only the brightest flares. 
However it seems unlikely that the flares observed at t > 10 5 



s are compatible with the extrapolation of the behaviour t 
found for early time flares; 

- the decoupling of the evolution of the peak luminosity of 
late time flares from the underlying continuum can account 
for the fainter distribution of the peak flux to continuum ratio 
AF/F of late time flares; 

- a sizable fraction of late time flares are compatible with vari- 
ability in the GRB afterglow. 

The global pictures that emerges is that the morphology of 
the flare light curve is the same for both early and late time 
flares, while they differ energetically. The similarities of late 
time flares with the early time ones and with the prompt emis- 
sion pulses as well as their energetic can be fairly well ex- 
plained in the framework of the accretion models. The exis- 
tence of flares i n this context is ascribed to instabilities ei- 
ther in the disk (|Perna et ai]|2006l) or in the fall-back material 
dKumar et alJ r2008; Rosswog 2007). This fragmentation can in 
principle account for the observed energetic and the longer du - 
ration of the accretion episodes at later times dPerna et al . 2006). 
Intermittence can be achieved also in presence of magnetic fields 
dProga & Begermanll2003t iProga & Zhangll2006h . Different ex- 
planations are provided by the magnetic reconne ction model 
dLvutikovll2006t lGianniosll2006t lLazaret all 120091) that do not 
require a long lasting activity of the central engine. However, the 
presence of ~ 86% of the late time flares of our sample that can 
be ascribed to variability in the GRB afterglow emission leads 
us to conclude that there could be no unique explanation about 
the nature of late time flares. 
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